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Abstract 

We calculate the contribution of the instanton - anti-instanton (II) pairs to 
the vacuum energy of QCD-like theories with Nf light fermions using the saddle 
point method. We find a qualitative change of the behavior: for Nf > 6 it starts to 
oscillate with Nf. Similar behaviour was known for quantum mechanical systems 
interacting with fermions. We discuss the possible consequences of this phenomenon, 
and its relation to the mechanism of chiral symmetry breaking in these theories. 
We also discuss the asymptotics of the perturbative series associated with the II 
contribution, comparing our results with those in literature. 
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In pure gauge theories the tunneling between topologically distinct classical 
vacua (described semiclassically by instantons) are known to shift the ground state 
energy down, as in quantum mechanics. It is often assumed that any non-perturbative 
effects should do the same: e.g. the bag constant of the MIT bag model was always 
assumed to be positive, without much discussion. However this is no longer obvi- 
ous if there are fermions in the theory. Moreover, in supersymmetric theories only 
positive shifts, if any, are allowed. 

In order to see how transition from one regime to another may happen, Balitsky 
and Yung |Q] considered a toy model, a quantum mechanical double-well potential 
coupled to a fermion field with a certain coupling constant a. The contribution of 
tunneling-antitunneling paths was found to have a complex phase proportional to 
a, which exactly flips sign at the a values at which the model becomes supersym- 
metricP]. 

In this work we study what happens in the QCD-like gauge theories with Nf 
light quarks. Unlike in the toy model mentioned, in this case the instanton con- 
tribution is only a part of the nonperturbative contribution to the vacuum energy. 
However, there are multiple evidences (see e.g. recent review 0]) that at least for 
Nf = — 3 (here and below we imply that the number of colors is N c = 3) it is 
a very important if not dominant part. The so called "instanton liquid" model || 
picture the vacuum of those theories as a relatively dilute and uncorrelated ensemble 
of instantons: it reproduces many phenomenological facts about hadronic correla- 
tions functions and spectroscopy, was recently directly supported by multiple lattice 
studies, see e.g. Q]. With increasing temperature T correlations between instantons 
and anti- instantons become more important, so that the random picture is no longer 
valid. Eventually a chiral restoration phase transition takes place, in which only 

strongly correlated instanton - anti-instanton pairs, or "molecules" || are present. 

The absolute value of also coincides with the known shift, which ensure that the method is 
correct. 
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This idea reproduces lattice data for close-to-QCD theories (Nf = — 3) on the 
position and type of the chiral phase transitions, hadronic screening masses etc. 

The pattern and mechanism of the chiral symmetry breaking at larger Nf is 
not yet understood. It is known that chiral symmetry should be restored in the 
vacuum above certain critical value Nf, which should be below the value where the 
asymptotic freedom disappears Nf < N^ F = l\N c /2. Between them the theory is 
conformal due to the infrared fixed point || . Studies of the interacting instanton en- 
semble H have shown that instantons alone cannot support the condensate already 



for Nf > 4. Lattice simulations for Nf = 8 [15] has found at weak coupling regime 
a chirally symmetric phase, similar to what was recently observed for Nf = 16 |14}] . 
Furthermore, for Nf = 4 [10] have the quark condensate was found to be nonzero 
but drastically smaller than for Nf = — 3. However, estimates based on one loop 
gap equation Jjj] suggest much larger critical value JVS ~ 11: presumably the con- 
densate induced by this mechanism is too small be seen in the lattice measurements 
mentioned. 

As Nf is increased, the two basic components of the instanton ensemble men- 
tioned above, (i) "single" (uncorrelated) ones and (ii) strongly correlated II pairs 
or "molecules", have the opposite trends. Singles can only exist due to nonzero 
value of the quark condensate. The fermionic determinant proportional to (| < 
qq > \p 3 ) Nf ) suppresses small-size instantons, so that only large ones, with sizes 
P ~ I < QQ > l ^ 1 ^ 3 ) survive^. 

The "molecules" do not need a nonzero quark condensate: thus they become the 
dominant effect. Furthermore, as noticed in |llj] , for increasing Nf they move toward 
smaller sizes, and for Nf > lliV c /2 — 3 their density even becomes ultraviolet diver- 
gent0. The differential contribution to the normalized molecular partition function 

2 Those are outside the semi-classical domain, unless we are close to chiral restoration boundary 

§■ 

3 It can be seen on purely dimensional ground: the density of molecules is dn m ~ dpA 2b p 2b 5 , 
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of instantons and anti-instantons with small (in order to use semi-classical formulae) 
radii pi, pj « 1/A: 



d2 ^L = y(4) C 2 4 f ^ {R)d n(-p 2 \T ir (R,n)f) Nf exp(-S int (R,m, 
oi dp j J 



dpi 

where R is the separation in units of p and Q the relative orientation of the II pair. 
C is standard single instanton density: 

4.6exp(-1.86iV c ) 1 2Nc 

c = ^(N c -mN c -2)i 7 (So) exp( " 5o) ' (2) 

Up to two loops the relation between the instanton action Sq and size is 

Sr. 2b bl 

p = e -ir\-S + l\W, (3) 

where b = ^N c - |iV/, &i = ^iV^ - ^NcNf + ^. The so called overlap matrix 
element of the Dirac operator enteing here is 

T n = J dtd 3 x(j)\(x - zi, O/) p(f)j(x - zj, n T ) (4) 

with <p being fermionic zero modes of the insatnton and anti-instanton. 

The gauge interaction Si n t(R,Q)) in general depends on the gauge configura- 
tions used. For II configurations one cannot use the equations of motion because 
there is no notrivial minimum: one should use the Streamline equation instead 
0, |l2ft . However, we will only be interested in the saddle points at large enough 
R, so that well known dipole formula would be sufficient. 

The integrand in (1) has a maximum at R — > 0, so the integral appears to be 
dominated by the weak field configurations, which belong to perturbative sector. 
The way to separate the non-perturbative physics is to use the saddle point method, 
moving the contour in the complex R plane Q. Let us split the contour into 2 parts, 
one going from R = to R = too, corresponding to the perturbative contribution, 
which becomes UV divergent at b < 2. This phenomenon is similar to the UV divergence in the 
0(3) a model: both are curious examples of the UV divergences of nonperturbative nature. 
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and the other one going from that point to R = +00 in such a way that \R\ is 
always big, so that (4) is applicable?]. The second integral is identified as a non- 
perturbative part of Z mo i. It is generally complex, and its imaginary part should 
be compensated by the perturbative contribution: thus the well known relation to 
perturbative series: 

E pert = _\ J°° AE_j m{E non. P ert {g)l (g) 

At the saddle point there is a balance between the gauge and fermionic part of 
the action: therefore we treat Nf as a largeparameter, which however should also be 
much smaller than the single instanton action Sq in order to ensure that the saddle 
point is at large R: So » Nf >> 1. If so, we may use asymptotic expressions for 
Tjj and Si n t fj Changing to variables (3 = l/R 2 and £ = cos 2 (6), where 9 is the 
only relevant orientation angle, (4) becomes: 

j2 ynon.pert 

a A mol _ (j2p6j 



I 



dhxpidhxpjVW 

16 cos 2 (6) 



Jo vr/2 J \ 



R e 



4(-W) N f £ dJ | o °° ^ exp(-5 4/3 2 (l - 4£) + (3JV> - 2) In/?). (6) 



There are two saddle points at /3q = ±y 8ffo(i-4g) • When £ < 1/4 they are real and 
we take the positive one, but for £ > 1/4 they are imaginary. It is irrelevant which 
one we choose (the choice is related to the definition of the perturbation series), so 
we take the one in the lower complex plane. Doing the (3 integral by the steepest 



Note that the addition of more terms in the expansion of Tjj and Sint leads to new saddle 



4 Presumably there are no singularities that can prevent it 

of 1 U 

points at smaller separations. Those however are artefacts of the truncation of the asymptotic 
series and should be disregarded. 
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descent method (the answer is analytic in £), we get for (8): 



3iV f -2 3N f -3 n M _ t A7 3iV f -2 

I = A(-16) N f^(8S e) &— (3JV>-2)^— ^ ^W^^(l-4£) ~ (7) 

The above integral has a singularity at £ = 1/4^| which can be avoided by taking 
£ integration over a contour in the complex plane. A contour in the upper £ plane 
matches the real positive /3q for £ < 1/4 to the negative imaginary one we have 
chosen for £ > 1/4. 

For even Nf the integral (10) can be written as a contour integral around the cut 
from to 1 (avoiding £ = 1/4), and the contributions are from two poles: £ = 1/4, 
which is of order (3/2)Nf — 1 and gives purely imaginary contribution and £ = oo, 
which is of order 3 — 1/2JV/, which gives a real contribution, but ceases to exist 
for Nf > 6. That is where the behaviour of the integral changes. For odd Nf one 
can express the integral via the incomplete elliptic integrals and their derivatives. 



In this case there are both real and imaginary contributions. Table 1. contains 
the values 



d 2 E 

the values of the real and imaginary parts of ^HEp/fln^ in units A, for N c = 3, for 



Table 1. 



6 It is actually an artefact, due to the truncation of (3 = 1/R 2 series. In fact close to that 
point the steepest descent method is not applicable if we truncate the gauge action to order 0(/? 3 ), 
because there is no large parameter in the exponent. However the 0(/3 4 ) term is nonzero at £ = 1/4 
and if we take it into account, the integral will be finite. In this case, however, we have 4 saddle 

points, and we have the difficulty of defining the integral so that no spurious ones contribute to it. 

7 Note that for small Nf the results are unreliable because the condition Nf » 1 is not true; 

those are shown for comparison only. Note also that we have not attempted to integrate over the 
instanton sizes, because that cannot be done without introducing a particular assumptions about 
the mechanism of the cutoff at large sizes. 
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Nf 


R,C ^ ^mol.gas 

dlnpjdlnpy 


j2 TT1 

Itti mol -s as 

' din pidln pj 


I(Nf) 


1 


.7959 * 10~ 5 


-.1717* 10- 4 


-.7173 * 10~ 5 


2 


0. 


-.5098 * 10~ 6 


-.2690 * 10~ 5 


3 


-.3668* 10- 7 


-.5474 * 10~ 7 


-.2833 * 10~ 5 


4 


-.1917* 10- 7 


-.5330 * 10- 8 


-.1991 * 10~ 5 


5 


-.1027* 10" 7 


.6464 * 10~ 8 


.1178* 10' 4 


6 


0. 


.1520* 10" 7 


.7913 * 10~ 4 


7 


.5067 * 10~ 7 


.1127* 10" 7 


.7147* 10" 4 


8 


0. 


-.2714* 10" 5 


-.2292 * 10~ 2 


9 


-.1898 * 10~ 3 


.1368* 10^ 4 


.2067 * 10~ 2 


10 


0. 


.3365 * 10- 4 


.4120 


11 


.1093 * 10~ 6 


.2421 * 10~ 8 


.1480 


12 


0. 


-.6030* 10" 11 


-132.8 



Here I(N f ) = Im f^J^ e 2S ° S 4jVc+(3/2)A ^ 1 i s the imaginary part of the 
energy derivative with the dependence on So (respectively g) removed. 

The obtained Nf dependence of the instanton- anti-instanton contribution to 
the partition function is very peculiar. First of all, in many cases the real part 
of the obtained result vanishes. Second, its signs start to oscillate for Nf > 5. 
Coincidentally, this is also the region where the instantons cannot support chiral 
condensate, and, as discussed above, their contribution is probably very small. 

In this paper we have not investigated the effective interaction between fermions 
induced by the "molecules" , but comment that one should expect for intermediate 
Nf appearance of non-perturbative forces and hadronic states having two distinct 
scales, dictated by the condensate and "molecules". Furthermore (at least in the 
saddle point approximation used in this work) their sign and magnitude of the 
latter should be proportional to the vacuum shifts evaluated above. As a result, 
the oscillatory behaviour found in this work should also propagate to correlators (at 
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sufficient small distances) and hadronic spectra. 

Let us now evaluate the large order coefficients in the perturbation series due 
instanton- anti-instanton contribution. Using (7) we can find the energy: 



Zvpert 1 d 87T 2*N c -(3/2)Nf + l _^ 

P =-- TfcTT^ e ^ T ( N f) 



d l E^ L 4= _1 r dg 8vr 
d\n pid\\i vr J g k+l g 2 

= ~m)± (s^)* 72 jf (S r^ 3/2)N ? +k/2 e- 2S °, (8) 

This is yet another saddle point integral with a saddle point at Sq = 2N C — 3/4iVy + 
k/4. The applicability condition now is k >> 3Nf. Finally we get the expected 
factorial behaviour^]: 

d2E T rt 4 T(AT \ k/2 



din ^id In p2 V 4-7T \87r 2 , 

4iV c -(3/2)iV / + fc /2 X ^-(3/2)iV /+fc / 2+ i/ 2 {3/2)Nf+k/2) 



2 



0) 



One can compare it with the results for the II contribution to cross section of 
the process of e + e~ — > hadrons. In the factorial behavior is (4iV c + fc/2)! for 



iVy = N c = 3. In [19] for the standard choices of iV c and Nf the authors get 



Re+e-^hadrons = £ -813(3280.5fc)- 35 / fc (10 + Wj^, (10) 

which is quite close to the behaviour of (11) for N c = 3 and Nf = 2. 

In summary, we have calculated the contribution of a correlated instanton- 
anti-instanton pair to the partition function of the QCD-like theories with increas- 
ing number of fermions. On general grounds (and the analogy to the quantum- 
mechanical problem with fermions) it was expected that the behaviour should have 

an oscillatory pattern, shifting the ground state down or upward as Nf grows. This 
8 The coefficients are smaller than those due to "renormalons" , which however correspond to a 
completely different set of diagrams, with maximal number of loops instead of classical "trees". 
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was indeed found to be true, but only for large enough Nf. The same is expected to 
happen for the hadronic spectra. Our results for the asymptotic of the perturbative 



series generated by II configurations are very close to those obtained in [18, by 
a different method. 
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